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ABSTRACT 

The distribution of the apparent total energy emitted by a gamma-ray burst 
reflects not only the distribution of the energy actually released by the burst 
engine, but also the distribution of beaming angles. Using the observed energy 
fluences, the detection thresholds and burst redshifts for three burst samples, I 
calculate the best-fit parameters for lognormal and power-law distributions of 
the apparent total energy. Two of the samples include a small number of bursts 
with spectroscopic redshifts, while the third sample has 220 bursts with redshifts 
determined by the proposed variability-luminosity correlation. I find different sets 
of parameter values for the three burst samples. The Bayesian odds ratio cannot 
distinguish between the two model distribution functions for the two smaller burst 
samples with spectroscopic redshifts, but does favor the lognormal distribution 
for the larger sample with variability-derived redshifts. The data do not rule 
out a distribution with a low energy tail which is currently unobservable. I find 
that neglecting the burst detection threshold biases the fitted distribution to be 
narrower with a higher average value than the true distribution; this demonstrates 
the importance of determining and reporting the effective detection threshold for 
bursts in a sample. 

Subject headings: gamma-rays: bursts 



1. Introduction 

The growing number of gamma-ray bursts with redshifts has not only established that 
most, if not all, bursts are at cosmological distances, and that up to 10 54 erg are radiated 
by a burst, but permits us to determine important intrinsic physical distributions. Since we 
sample a burst's radiation pattern at only one point, the observed energy fluence can only 
be related to the energy flux emitted in our direction; this energy flux can be expressed as 
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the total energy the burst would have emitted if it radiated isotropically. Here I consider 
the distribution of this apparent total energy. 

The distribution of the apparent total energy is a convolution of the distribution of the 
actual energy emitted and the distribution of the angle into which the emission is beamed, 
both quantities of crucial importance in understanding the physics of the progenitor, the 
frequency of burst occurrence, and the impact of a burst on its environment. Frail et al. 
(2001) have recently determined the beaming angles for a burst sample by modeling the 
late-time breaks in the temporal decay of the afterglows; accounting for the beaming angle 
provides the actual total energy, which they found is clustered around E ~ 5 x 10 50 erg. If 
the analysis of Frail et al. is indeed correct (there are competing models for the evolution of 
the afterglow) then methods similar to those I develop here will be necessary to determine 
properly the beaming angle and total energy distributions. 

This study extends the work of Jimenez, Band & Piran (2001) which considered lognor- 
mal distributions for the apparent total gamma-ray energy, the peak gamma-ray luminosity, 
and the total X-ray afterglow energy. Jimenez et al. expanded the database of bursts with 
spectroscopic redshifts by adding bursts for which a redshift probability distribution could 
be inferred from the host galaxy brightness. Currently the bursts with only redshift proba- 
bility distributions do not augment the relevant burst database sufficiently to warrant their 
inclusion in my study, although in the future it may be advantageous to include these bursts 
if the determination of spectroscopic redshifts does not keep pace with the detection of bursts 
and their host galaxies. 

Recently correlations have been proposed between burst properties and their peak lumi- 
nosities. Norris, Marani & Bonnel (1999) proposed that more luminous bursts have smaller 
time lags between energy channels, while Fenimore & Ramirez-Ruiz (2001) reported that 
the light curves of more luminous bursts were more variable; Schaefer, Deng & Band (2001) 
found that applying both correlations to the same burst sample resulted in consistent burst 
luminosities. The redshift is determined from the derived luminosities and the observed peak 
fluxes. Thus these correlations can give us a large burst sample with redshifts from which 
the energy distribution can be determined. 

In §2 I discuss the methodology for finding the best parameter values for a given func- 
tional form of the distribution, and for comparing different functional forms. I also evaluate 
the sensitivity to small burst samples, and demonstrate the importance of considering the 
detection threshold in determining the distributions. The energy distributions resulting from 
two small burst samples with spectroscopic redshifts and from a large sample with redshifts 
from the variability-luminosity correlation are described in §3. Finally, §4 summarizes my 
conclusions. 
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2. Methodology 

2.1. The Likelihood Function 

I begin with an assumed distribution p(E | dj, Mj, I) where E is the apparent total burst 
energy, d} is the set of parameters which characterize the jth model distribution function 
represented by Mj, and / specifies general assumptions about the distribution function. I use 
p(a | b) to mean the probability of a given b. Below I present the distribution functions used 
in this study. A distribution is assumed to be universal and not a function of redshift (i.e., 
no evolution), or of properties of the host galaxy, burst, etc. Of course, once a sufficiently 
large sample is available, the energy distributions for burst subsets can be investigated. 
I use normalized distributions since the burst rate is not of interest here. The observed 
energy fluence F and the burst energy E are related by F = E(l + z) / An D L (z) 2 = EC(z) 
where Dl(z) is the luminosity distance. It is in calculating the luminosity distance that the 
burst redshift and a cosmological model are required. Here I assume H =65 km s -1 Mpc" 1 , 
Vt m = 0.3 and fi A = 0.7. 

This energy distribution is converted into p(F \ dj, z, Mj, I), the probability of obtaining 
the energy fluence F given the parameters d} for model Mj, the burst redshift z, and other 
assumptions / (e.g., the choice of cosmological model). However, the observed fluences are 
not drawn from this probability distribution but from p(F | Ft, dj, z, Mj, I), the normalized 
distribution which is truncated below Ft, the minimum fluence at which that particular 
burst would have been included in the burst sample. 

I now form the likelihood 

N 

L j = Y[p(Fi I *T,i, %, Zi, Mj, /) (1) 

1=1 

where the ith burst has energy fluence Fi, threshold fluence Ft^, and redshift Z{. In the 
"frequentist" framework best fit parameters are typically found by maximizing Lj. 

The likelihood is Lj = p(D | d}, Mj, I), where D is the set of observed fluences and 
fluence thresholds. The Bayesian method of estimating the parameters is to calculate 

(a-> = j d a,a lP[ S,\D,M,J) ^ (2) 

/ ddj <Tj p(D | dj, Mj, I)p(dj \ Mj, I) _ J da} d* Yl i=1 p(Fi \ F Tji , d}, z i} Mj, I)p{aj \ Mj, I) 
$ da j p{D\a j ,M j ,I)p{a j \M j ,I) ~ J da 3 Y^ =l p{F l \F T , l ,<i ] ,z l ,M 3 ,I)p{ ( x 3 \M 1 ,I) 

where p(d} | Mj, I) is the prior for the parameters d}. If A,- = p(D \ dj, Mj, I)p(dj \ Mj, I) is 
sharply peaked, then the expectation value of the parameters occurs at the peak of A^. Note 
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that Aj is the likelihood in eq. 1 times the priors for the parameters. In some sense the 
priors indicate the natural variables in terms of which the priors are constant. For example, 
here the burst energies vary over a number of decades, and thus I assume that the priors 
for average or cutoff energies are logarithmic. Consequently maximizing Aj is equivalent to 
maximizing the likelihood in terms of the logarithm of the energy parameters. This is the 
methodology I use here; I do not attempt to integrate the integrals in eq. 2. 

For each burst the cumulative probability is 



If the assumed energy distribution function is an acceptable characterization of the ob- 
servations (which would be the case if the model Mj is correct) and all the assumptions 
are valid (e.g., the cosmological model is correct), then the cumulative probabilities -P(Fj) 
should be uniformly distributed between and 1, and have an average value of (P(F i )) = 
1/2 ± (12iV)~ 1/2 for N bursts in the sample. 



As mentioned above, the average of the cumulative probability (P(Fj)) should be 1/2 
(within a quantifiable uncertainty) if the distribution function describes the observations 
satisfactorily. Thus the value of this statistic for model Mj is a measure of the acceptability 
of that model. Further, a comparison of the values for different models is a measure of the 
relative merits of these models; of course, this comparison should account for the expected 
uncertainty in the value of this statistic. 

The Bayesian framework provides a clear prescription for comparing models through the 
odds ratio. Let p(Mj | D, I) be the posterior probability for the jth model Mj given the data 
D. Then the odds ratio comparing the jth and kth models is Ojk = p{Mj | D, I)/p(Mk \D,I). 
But by Bayes' Theorem p(Mj \ D, I) oc p(Mj \ I)p(D \ Mj, I) where the normalizing fac- 
tor is independent of Mj, and thus cancels in forming the odds ratio. I assume that no 
model is favored a priori; therefore the "priors" p(Mj \ I) are the same for all Mj, and 
cancel in forming the odds ratio. Here D consists of the observed energy fluences (and 
the detection thresholds). Thus p(Mj \ D, I) oc YliPi^i] F T ,i, z i: Mj, I). However, we be- 
gin with the more fundamental probabilities p(Fi\ Fx,i, Zi, Mj,ctj, I), which are functions 
of the model parameters dj. In determining the preferred model we are uninterested in 
the specific model parameters, and therefore we "marginalize" over the parameters d}: 



p(Fi | F T ,i, z u Mj, I) = / d<Tjp(aj | Mj, I) fliP^ | F T>i , z h Mj, d}, /), where p(ctj \ Mj, I) is the 




(3) 



2.2. Model Comparison 
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prior for the parameters of the jih model. Thus 

, k = I d6: i Qf P(Fi I Fr,i, Zj, Mj, a}, i>(d} \ M,, I) 
1 J da k f]f p(Fi | F Tji , Zi, Mi... a k , I)p{a k \ M k , I) 

It will be noted that the odds ratio is the ratio of the likelihoods for each model (eq. 1) 
marginalized over the model parameters. Often the prior for the parameters is set equal to 
a delta function at the parameter values which maximize the likelihood, which is of course 
circular reasoning. In this common formulation, the odds ratio is the ratio of the peak values 
of the likelihoods for each model. Below I present priors defined in terms of the expected 
parameter ranges, in which case the integrals in eq. 4 must be integrated. I present values of 
the odds ratio with both the delta function priors and the more refined priors. It is possible 
that with the delta function priors the odds ratio may favor one model, but with the more 
refined priors another model is favored, as is indeed the case here. 



2.3. Distribution Functions 



2.3.1. Lognormal Distribution 

I assume that E has a log-normal distribution 

(\nE -\nE) 2 



p(E\E ,a)d(\nE) = 



2na 



cxp 



2a 2 



d(lnE) . 
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Thus the fluence F also has a log-normal distribution. Note that a is a width in logarith- 
mic space, and the linear change of variables from E to F does not affect this width. As 
discussed above, we need to consider the range over which the fluence could actually have 
been observed, i.e., for fluences above the threshold F T . The resulting normalized fluence 
probability distribution is 



p ohs (F\F T ,E ,a,z)d(\nF) 





(\n[E C(z)]-lnF) 2 ' 


0(F- 


F T )d(\nF) 




[ 2a2 


\ (l + erf 


"ln[£ C(z)]-ln(F T )" 
V2a 


) 



(6) 



where Q(x) is the Heaviside function (1 above x — 0, and below) and C(z) = (1 + 
z) I '4-7T D l{z) 2 converts energies in the burst's frame to fluences in our frame; the denominator 
results from integrating over the numerator from Ft to infinity. 



For the Bayesian formulation we also need the priors for the model parameters. There 
is no reason to favor one energy over another over many energy decades, and thus I assume the 
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prior is constant in logarithimic space: p(E )dE = (E Yn[E u / E^^dE^ = (\n[E u / Ei])~ x d In E = 
{\og 1Q [E u / Ei])~ l d\og l0 EQ where E u and E t are the upper and lower limits of the permitted 
range (because of the logarithmic dependence, the result is not very sensitive to the precise 
values). I will use E\ = 10 51 erg and E u = 10 54 erg. Similarly, I have no a priori information 
about a, and therefore assign it a uniform prior between and 5. Note that the distribution 
function is explicitly a function of In E$ and a, and the priors indicate that these are indeed 
the natural parameters. 



2.3.2. Single Component Power Law Distribution 

I assume that the underlying energy probability distribution is 

p(E\E 1 ,E 2 ,a)dE = J^S_^E- a dE ; E 1 < E < E 2 , a^l 

T^m dE 5 E,<E<E 2 , a = l . (7) 

If Ei is extended to or E 2 to infinity, then we must restrict a to be less than or greater 
than 1, respectively. The expected fluence probability distribution, accounting for the fluence 
threshold F T) is 

p(F\F T , El ,E 2 ,a)dF= ^K'gfi F-°dF ; E 1 < F/C(z) < E 2 , a^l 

where again C(z) = (1 + z) / An D L (z) 2 converts burst energies to fluences. 

As with E for the lognormal distribution, I have no reason to prefer any value of the 
energy limits E 1 and E 2 over many energy decades, and therefore I again use logarithmic 
priors for these two energies. For definiteness, I will assume that E 1 can have a value between 
10 49 and 10 52 erg, and E 2 between 10 52 and 10 55 erg. The spectral index is assumed to have 
a uniform prior between —2.5 and 2.5. 



2.4. Data 

The methodology discussed above requires the energy fluence F, the threshold fluence 
F T and the redshift z for each burst. I consider 3 samples. The B9 sample consists of 9 
bursts with BATSE data and spectroscopic redshifts. The fluences were calculated by fitting 
the BATSE spectrum accumulated over the entire burst with the "GRB" function (Band 
et al. 1993), and then integrating the resulting fit over the 20-2000 keV energy range in 
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the burst's rest frame and over the time during which the spectrum was accumulated. The 
resulting fits are presented in Jimenez et al. (2001). There are some bursts for which the 
high energy power law in the GRB function has an index (3 < —2 (where N oc E 13 ) and thus 
the ffuence depends crucially on the high energy cutoff (or roll-off) which must exist for a 
finite ffuence but which could not be determined from the BATSE data. The spectroscopic 
redshifts are taken from Frail et al. (2001). 

The limiting ffuence Ft is more difficult to determine. The bursts in our sample must 
have been intense enough to be first detected and then properly localized. In addition, a 
decision was made to attempt to observe the afterglow and thus determine the redshift. These 
threshold quantities have generally not been reported. Note also that detectors almost never 
trigger on the ffuence but usually trigger on the peak count rate sampled over a time bin of 
a specified duration. Here I will assume that the ratio of the observed to threshold ffuence, 
F/Ft, is the same as the ratio of the observed to threshold peak count rate, C max /C min , for 
the BATSE data. These thresholds are most likely underestimates of the true thresholds. 
The C max /C m i n ratio is a standard part of the BATSE catalog (BATSE team 2001); however, 
in some cases the BATSE team did not calculate this quantity because of data gaps, in which 
case 1 estimated this ratio from the lightcurves. This sample is described by Table 1. 

The C17 sample consists of the f 7 bursts with spectroscopic redshifts and fluences in 
Frail et al. (200f). This sample is basically a superset of the B9 sample with an addi- 
tional BATSE burst for which there are no spectra, and bursts observed by Beppo-SAX and 
Ulysses. For those bursts without reported detection thresholds I use a ffuence threshold of 
f0~ 6 erg cm~ 2 . 

Finally, the F220 sample uses the 220 bursts in Fenimore & Ramirez-Ruiz (200f ) with 
redshifts determined by the variability-luminosity correlation. This sample was selected 
to have a peak count rate accumulated on the 256 millisecond timescale of greater than 
f.5 cts s _1 , which provides a well-defined detection threshold. The fluences were taken from 
the BATSE catalog without any k- correct ions. Note that Bloom et al. (200f ) find that the 
k-correction for the 20-2000 keV energy range is of order unity. 

It should be noted that for the first two samples the redshifts are reliable but the 
detection threshold is very uncertain. Even when the detection threshold is known for the 
gamma-ray portion of the burst, the effective threshold for optical follow-up observations 
has not been reported. On the other hand, the detection threshold for the third sample is 
known, but the validity of the variability-luminosity correlation is still not well established, 
and the uncertainty in the resulting redshifts is not considered. 
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2.5. Simulations 

To determine the sensitivity of the methodology to the number of bursts and to demon- 
strate the importance of considering the fluence threshold, I ran a series of simulations. For 
each simulation I first created between 100 and 500 simulated databases to which I then 
applied the methodology described above to determine the parameters of their energy dis- 
tribution. For some simulations I found the parameters with the fluence thresholds used in 
creating the database or with much smaller thresholds. 

For each simulated burst I needed a redshift, a burst energy and a fluence threshold. 
The redshifts were drawn from a distribution which is similar to the proposed cosmic star 
formation rate. The fluence threshold was drawn from a uniform logarithmic distribution 
over one decade. The burst energy was drawn from a lognormal distribution with specified 
central energy and logarithmic width, as long as the resulting fluence was greater than the 
fluence threshold. 

Figure 1 shows the importance of considering the fluence threshold in calculating the 
likelihood. As can be seen, the best fitted parameter values cluster around the input central 
energy Eq = 10 53 erg and logarithmic width a = 1.5 when the fluence threshold is considered 
(asterisks), but cluster around a higher energy and narrower width when the fluence threshold 
is neglected (diamonds). The fluence threshold removes low energy bursts, resulting in a 
narrower apparent distribution which is shifted to higher energy. Each of the 100 simulated 
datasets had 80 bursts, and a fluence threshold between 10~ 6 and 10~ 5 erg cm -2 . Figure 2 
shows that the likelihood contours for a sample analyzed with the correct fluence thresholds 
(left panel) and thresholds a factor of 10 smaller (right panel) differ significantly. 

Table 2 gives the width of the distributions of the parameters of the lognormal distribu- 
tion, logi?o and a, for databases with 9, 20, 40 and 80 bursts. As expected, the distributions 
become narrower as the number of bursts increases. As can be seen, a database with 40 
bursts should give satisfactory best-fit parameter values. 

3. Results 

As can be seen from Table 3, the parameter values at the likelihood maximum for the 
two distribution functions differ for the three burst samples, and the 90% confidence ranges 
from one burst sample do not always include the parameters from the other samples. Note 
that the lower energy cutoff E x for the simple power law distribution is not fitted. The 
fits are insensitive to E 1 values smaller than the lowest energy for which any burst in the 
sample would have been detected, and E 1 cannot be greater than the smallest observed 



- 9- 



burst energy; the difference between these two limits is very small. Figures 3-8 show the 
likelihood contours for the fitted parameters. A ridge of somewhat lower likelihood values 
curves towards lower values of E and higher values of a for the lognormal distributions. 
Thus the data do not strongly exclude a broader distribution which includes lower energy 
bursts which are not detected because their fluence is below the threshold. Consequently, the 
90% confidence bounds for the parameters are quite broad. For the power law distributions, 
the favored high energy cutoff is at the highest observed energy in the burst sample; if the 
bursts truely have a power law distribution, this cutoff is most likely somewhat greater. 

Table 3 also provides the average of the cumulative probabilities for each sample and 
distribution. As can be seen, these averages are within ~ la of the expected value of 
1/2, indicating that the distributions are acceptable descriptions of the data, and that this 
statistic cannot be used to discriminate between distributions. The actual distributions of 
these cumulative probabilities are also presented by Figures 3-8. 

Finally, Table 3 presents the likelihood ratios and the Bayesian odds ratios comparing 
the lognormal to power law distributions. The likelihood ratio is a Bayesian odds ratio using 
delta function priors set to the parameters which maximize the likelihoods (which violates 
the definition of a prior). The Bayesian odds ratios given here use the priors described in 
§2.3. The ratios show that using a broad prior distribution favors the lognormal distribution 
over the power law distribution. Based on the odds ratio, the B9 and C17 samples are 
insufficient to discriminate between the two distribution functions. On the other hand, the 
odds ratio favors the lognormal distribution for the F220 sample. 

4. Discussion and Conclusions 

The parameters and parameter ranges differ for the three different burst samples. It is 
not clear whether we yet have a sufficiently large, properly defined burst sample from which to 
calculate the energy distribution. The two samples with spectroscopic redshifts do not have 
correct detection thresholds: the threshold for detecting the burst itself is usually reported, 
but the intensity threshold which triggers further localization and spectroscopic redshift 
determination has not been reported. Indeed, there may not yet be a formal definition of such 
a followup threshold. Thus these two samples are flawed. On the other hand, the validity of 
the variability-determined redshifts has not yet been proven, although the detection threshold 
was defined in choosing the burst sample. The importance of the detection thresholds for 
statistical studies of burst samples argues for well-defined (and reported) intensity thresholds 
for triggering the followup observations of the expected large number of HETE-II and SWIFT 
burst localizations. 
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The distributions of cumulative probabilities and the averages of these distributions 
indicate that the two functional forms used here are sufficient to describe the distribution 
of energies; consequently, I concluded that the data do not justify trying more complicated 
distribution functions at this time. The Bayesian odds ratio does not distinguish between 
these two functional forms for the two samples with spectroscopic redshifts, although it does 
favor the lognormal distribution for the large F220 sample with redshifts derived from the 
variability of the burst lightcurves. 

The energy distribution cannot be determined for energies below the lowest energy 
threshold (i.e., the lowest burst energy corresponding to the ffuence thresholds). Indeed, 
Hakkila et al. (1996) make a distinction between the "observed" and "intrinsic" luminosity 
functions in their study of luminosity functions for cosmological bursts based on the shape 
of the peak flux distribution; they point out that the observed distribution may be much 
narrower than the intrinsic distribution. In my study the true low energy cutoff of the power 
law distribution cannot be determined. Similarly, the likelihood contours for the lognormal 
distribution do not rule out broader distributions with lower central energies. 

We anticipate that the HETE-II and S WIFT missions will result in the construction of 
a large burst sample with spectroscopic redshifts. Properly defined subsets can be studied 
to identify trends with burst redshift, duration, and other properties. As we move from the 
study of individual bursts to the study of burst ensembles, we must define and report the 
criteria (e.g., detection thresholds) by which the burst samples are collected. 

This work was performed under the auspices of the U.S. Department of Energy by the 
Los Alamos National Laboratory under Contract No. W-7405-Eng-36. 
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Fig. 1. — Best-fit parameter values including (asterisks) or neglecting (diamonds) the fluence 
threshold in calculating the likelihood for 100 simulated databases with 80 bursts each. The 
model lognormal energy distribution had a central energy of E = 10 53 erg and a logarithmic 
width of a = 1.5. The median values for the fits including (large solid cross) or neglecting 
(large dashed cross) the fluence threshold are indicated. The bursts were drawn from a 
redshift distribution similar to that of star formation, and the fluence threshold was between 
erg cm" 



10~ 6 and 10~ 5 — — ~ 2 
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Fig. 2. — Comparison of the likelihood contours for a sample analyzed with the correct 
ffuence thresholds (left) and thresholds a factor of 10 smaller (right). The sample of 9 
bursts was drawn from a lognormal distribution with Eq = 10 53 erg and a = 1.5, a redshift 
distribution similar to star formation, and a uniform fluence threshold between 1CT 5 and 
10~ 4 erg cm -2 . 
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Fig. 3. — Contour plot of the likelihood for the lognormal energy distribution for the B9 
sample (left panel). The parameters are the central energy Eq and the logarithmic width a. 
The asterisk indicates the location of the maximum likelihood while contours are spaced by 
Alog 10 L = 0.1 starting from the maximum value. Cumulative distribution of the cumula- 
tive probability for each burst assuming their energies are drawn from the best-fit lognormal 
energy distribution for the B9 sample (right panel). The dashed line is the expected distri- 
bution. 
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Fig. 4. — The same as Figure 3, but for a simple power law energy distribution for the B9 
sample. The parameters are the upper cutoff energy E 2 and the power law index a. 
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Fig. 5. — The same as Figure 3 for the lognormal distribution and the C17 sample. 
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Fig. 6. — The same as Figure 4 for the power law distribution and the C17 sample. 
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Fig. 7. — The same as Figures 3 and 5 for the lognormal distribution and the F220 sample. 
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Fig. 8. — The same as Figures 4 and 6 for the power law distribution and the F220 sample. 
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Table 1. The BATSE Gamma-Ray Bursts Sample. 



Burst 


TP a 
r obs 

(erg cm _/ 




TP b 
r burst 

(erg cm~ z 




^obs 


Cm ax 
C m i n 


-Eobs° 

(10 51 erg) 


TP d 
-^burst 

(10 51 erg) 


970508 


3.18 x 10~ 


-6 


2.59 x 10- 


-6 


0.835 


3.3 f 


6.734 


5.482 


970828 


9.57 x 10' 


-5 


7.88 x 10- 


-5 


0.958 


20 f 


267.1 


219.3 


971214 


9.44 x 10- 


-6 


7.59 x 10- 


-6 


3.412 


7.32° 


261.3 


210.7 


980703 


2.26 x 10- 


-5 


2.13 x 10- 


-5 


0.966 


3.08° 


63.94 


60.18 


990123 


2.68 x 10- 


-4 


1.93 x 10- 


-4 


1.600 


80.1° 


1996. 


1438. 


990506 


1.94 x 10- 


-4 


1.69 x 10- 


-4 


1.2 


50 f 


838.8 


854.0 


990510 


2.26 x 10- 


-5 


2.32 x 10- 


-5 


1.619 


19.3° 


172.0 


176.8 


991216 


1.93 x 10- 


-4 


1.70 x 10- 


-4 


1.02 


144. c 


611.1 


534.1 


000131 


4.18 x 10- 


-5 


2.71 x 10- 


-5 


4.5 


3 f 


1791. 


1159. 



a Fluence over 20-2000 keV in the observer's frame. 
b Fluence over 20-2000 keV in the burst's frame. 

c Gamma-ray energy over 20-2000 keV in the observer's frame, assuming 
isotropic emission, H = 65 km s^ 1 Mpc -1 , il M =0.3 and f2A=0.7. 

d Gamma-ray energy over 20-2000 keV in the burst's frame. 
e From the online BATSE catalog. 
f Estimated from lightcurve. 
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Table 2. Width of Parameter Distribution 



Number of Bursts 


Width a oflogE 


Width a of a 


in Sample 


Distribution 


Distribution 


9 


0.38 


0.47 


20 


0.40 


0.52 


40 


0.20 


0.31 


80 


0.13 


0.19 



Note. — In these simulations 100 samples were con- 
structed with the indicated number of bursts per sam- 
ple. The burst energies were drawn from a lognormal 
distribution with central energy Eq = 10 53 erg and loga- 
rithmic width a = 1.5. The redshift distribution mimics 
the cosmic star formation rate, and the threshold flu- 
ence was between 10 -6 and 10~ 5 erg cm -2 . The best fit 
parameters were found by maximizing the likelihood. 

a The width given is the range within which 1/2 the 
simulated bursts fell. 
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Table 3. Comparison of Model Distributions 



Quantity 


B9 a 


C17 b 


F220 c 




125.9 


52.48 


11.79 




1.6-315 


1.6-100 


2.-23.4 




1.87 


2.06 


1.88 




1.5-4.5 


1.7-4.15 


1.65-2.3 




1.6 


0.55 


0.12 




1440 


1460 


5000 




1440-3550 


1460-3350 


5000-6100 


a h 


0.74 


0.96 


1.81 




0.4-1.2 


0.75-1.25 


1.7-1.94 


(P)ln l 


0.4525 


0.4638 


0.4753 


(P)pl i 


0.4608 


0.4723 


0.4892 


r, k 

a (p) 


0.0962 


0.0700 


0.0195 


Likelihood Ratio 1 


4.29 x 10~ 2 


5.61 x 10~ 2 


4.38 x 10 2 


Odds Ratio 111 


1.19 


1.65 


9.92 x 10 3 



a Sample of 9 BATSE bursts with spectroscopic rcdshifts and fitted spectra (Table 1). 
b Sample of 17 bursts with spectroscopic redshifts (Frail et al. 2001). 

c Sample of 220 bursts with redshifts derived from variability redshifts (Fenimore & Ramirez- Ruiz 2001). 

d The central energy of the lognormal distribution, in units of 10 51 erg. The following line in the table 
provides the 90% confidence range. 

°The logarithmic width (in units of the energy's natural logarithm) for the lognormal distribution. The 
following line provides the 90% confidence range. 

f The low energy cutoff of the power law distribution, in units of 10 51 erg. This energy has been set to the 
lowest threshold energy for the sample. 

s The high energy cutoff of the power law distribution, in units of 10 51 erg. The following line provides the 
90% confidence range. 

h The power law index of the power law distribution, p(E) oc E~ a . The following line provides the 90% 
confidence range. 

'Average of the cumulative probabilities for the lognormal distribution; 1/2 is expected. 
j Average of the cumulative probabilities for the power law distribution; 1/2 is expected. 
k The standard deviation \/\\2N] 1 / 2 of the average of the cumulative probabilities for a sample of N bursts. 

'Ratio of the maximum likelihood for the lognormal to power law distributions. A value greater than 1 
favors the lognormal distribution. 

m Odds ratio comparing the lognormal to power law distributions. A value greater than 1 favors the 
lognormal distribution. 



